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Abstract 

We consider the evolution of an incompressible two-dimensional perfect fluid as the 
boundary of its domain is deformed in a prescribed fashion. The flow is taken to 
be initially steady, and the boundary deformation is assumed to be slow compared 
to the fluid motion. The Eulerian flow is found to remain approximately steady 
throughout the evolution. At leading order, the velocity field depends instanta- 
neously on the shape of the domain boundary, and it is determined by the steadi- 
ness and vorticity-preservation conditions. This is made explicit by reformulating 
the problem in terms of an area-preserving diffeomorphism g\ which transports the 
vorticity. The first-order correction to the velocity field is linear in the boundary 
velocity, and we show how it can be computed from the time-derivative of ^a- 

The evolution of the Lagrangian position of fluid particles is also examined. 
Thanks to vorticity conservation, this position can be specifled by an angle-like 
coordinate along vorticity contours. An evolution equation for this angle is derived, 
and the net change in angle resulting from a cyclic deformation of the domain 
boundary is calculated. This includes a geometric contribution which can be ex- 
pressed as the integral of a certain curvature over the interior of the circuit that is 
traced by the parameters defining the deforming boundary. 

A perturbation approach using Lie series is developed for the computation of 
both the Eulerian flow and geometric angle for small deformations of the boundary. 
Explicit results are presented for the evolution of nearly axisymmetric flows in 
slightly deformed discs. 
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1 Introduction 



This paper examines the dynamics of a two-dimensional (2d) fluid inside a 
container whose boundary is deformed slowly. The fluid is assumed to be per- 
fect and incompressible; consistent with the latter assumption, the area of 
the container is constant. Beyond potential applications such as the control 
of fluid flows, we use the problem as a paradigm for the study of Hamilto- 
nian fluid models depending on slowly varying parameters. This is an obvious 
flrst step: the 2d Euler e quations governing incompres sible perfect fluids are 



indeed Hamiltonian (e.g. . iMorrisoru . Il998l : ISalmonl . Il988l ) . and imposing bound- 



ary deformations is arguably the most natural way of introducing a parameter 
dependence. As is well-known in flnite dimensions, Hamiltonian systems are 
strongly constrained; as a result, slow changes of parameters lead to a re- 
mar kable behaviour encapsulated in the the ory of adiabatic invari ance (cf. , 



1985 



Arnoldl, Il989l: iLandau and Lifshitzl . Il960l ) and geometric angles (IHannayl . 
Berrvl. Il985f) . In 2d Euler, the material invariance of vorticity (e.g. 



Saffmanl . Il99i ) similarly imposes a strong contraint on the system, which we 



exploit extensively to derive what can be interpreted as fluid-dynamical ver- 
sions of adiabatic invariance and geometric angle. 



The problem we consider here is rather involved in its full generality. To make 
progress, we make a number of assumptions and consider the following sce- 
nario. At an initial time, a steady flow is given in some simply-connected 
domain Dq. The streamlines have the simplest topology, that of nested closed 
curves, and the flow is Arnold stable (see section 2 below). We then assume 
that this continues to hold throughout the evolution as the domain is being 
deformed. With these hypotheses, we use an asymptotic approach, based on 
the separation between the timescales of the boundary deformation and that 
of the flow, to answer two questions: (i) what is the leading-order approxi- 
mation to the (Eulerian) flow at any time; and (ii) what is the (Lagrangian) 
position of the fluid particles? 

The flrst question is answered by showing that the leading-order flow is steady 
at all times. This makes it possible to rephrase the problem in terms of an area- 
preserving diffeomorphism gj>^ which maps the vorticity in the initial domain 
to the vorticity in the deformed domain. The uniqueness of up to dis- 
placements along lines of constant vorticity, established in Wirosoetisno and 
Vanneste (2005; henceforth WV) and revisited here, shows that the leading- 
order velocity fleld is completely determined by the instantaneous shape of 
the boundary and is independent of the history of past shapes. This, of 
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course, is analogous to the adiabatic invariance of the action variables for 
finite-dimensional Hamiltonian systems with slowly varying parameters. 



To answer the second question, we need to compute the first-order correc- 
tion to the approximate velocity field obtained in (i). This is because the 
evolution equation for the particle position needs to be integrated over the 
long time scales required to achieve order-one boundary deformations. The 
first-order correction to the velocity field is linear in the boundary velocity, 
and it can be derived from dgf^/dt by solving a pseudodifferential equation. 
Once this is done, the particle-position problem reduces to the solution of (in- 
dependent) one-degree-of-freedom Hamiltonian systems with slowly varying 
parameters. Since particles remain on vorticity contours (which in effect are 
contours of constant action), only the position along each contour, regarded as 
an angle variable, needs to be determined. The value of this angle is found to 
depend on the history of the boundary shape. It includes a geometric contribu- 
tion, sim i lar to the Hannay-Berry angle, which possesses a nice interpretation 



Marsden et al 



(IHannavi. 11985 ; iBerryl . Il985l : iMontgomeryl . Il988l : IShapere and Wilczekl . Il989al : 



19901 ) ■ We note that the geon i etric angle h as been studied 



in fluid dynamics by IShashikanth and Newton! (119971 . 119991) who considered 
point- v ortex solutions of the 2d Euler equation, and by IShapere and Wilczek 
(jl989bl ) for Stokes flow; here it appears in the context of smooth inviscid flows. 



The determination of the leading-order Eulerian flow from the steadiness and 
vorticity-preservation conditions was treated in WV, where conditions for the 
existence of g\ and the uniqueness of the resulting velocity field were given 
in appropriate function spaces for sufficiently small boundary deformations. 
In the present article we adopt a more informal approach to treat both the 
Eulerian and Lagrangian problems under a slightly different set of hypotheses; 
rigorous proof of the adiabatic invariance will be the subject of a future work. 
It proves convenient to express our derivation in the language of differential 
forms in the space of the parameters defining the boundary shape. This makes 
explicit the linear dependence of several important quantities on the boundary 
velocity, and it gives a natural description of the geometric angle in terms of 
a curvature form in the parameter space. We use this language mainly as a 
notational tool, but it is clear that a more abstract geometric interpretation 
of the results could be given. This is discussed at the end of the paper. 

In the following section, we present a short description of the 2d Euler equa- 
tion in a deforming domain in order to fix the notation, and we consider the 
behaviour of the leading-order Eulerian flow. Next, in §3 we compute the 
first-order correction to the Eulerian flow, which depends only on the instan- 
taneous shape of the boundary and its velocity. Using these results, we study 
the Lagrangian flow in §4, where the geometric angle of the particle position is 
derived. In these sections, we consider general domains and arbitrary bound- 
ary deformations, requiring only that the boundary deformation be slow. The 
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Fig. 1. Parameterization of the shape of the domain Da by A: as the parameter A 
moves from A(0) = to A(ef) in £, the fluid's domain changes its shape from Do 
to D\, inducing a change in the leading-order Eulerian flow whose streamlines are 
indicated. 

results are given as solutions of partial (pseudo) differential equations, which 
in general will have to be solved numerically. In §5 we develop a perturbative 
approach for the solution of these equations, based on the assumption of small 
(total) boundary deformation. We carry out the calculation to second order, 
but the Lie series formulation that we employ is well suited for systematic 
extensions to higher orders. An application to nearly axisymmetric flows in a 
slightly deformed disc is presented in §6, followed by a Discussion. Technical 
details are relegated to the Appendices. 



2 Eulerian Flow: Adiabatic Invciriance 

We begin by studying the behaviour of the Eulerian flow. 
2. 1 Formulation 

Let -DA(£t) C be a simply-connected, bounded and smooth domain which is 
slowly evolving in time t in a prescribed fashion while keeping its area fixed. 
Here, A denotes the set of parameters defining the shape of the boundary, and 
the slowness of their time dependence is made explicit by the introduction 
of the asymptotic parameter < £ <^ 1. Denoting the (generally infinite- 
dimensional) space in which A lives by £, we can think of the evolution of 
the domain shape as the tracing of a curve h.{et) C C] see Figure 1. In this 
article, our concern is the behaviour of the flow as e ^ 0, that is, in the 
limit of slow boundary deformation, over 0{e^^) time scales so that 0(1) 
deformations are achieved. We make a blanket assumption that all functions 
are sufficiently smooth for our purposes, and we denote by C(Da) the space 
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of smooth real- valued functions in Da- 

We can describe the evolution of a perfect fluid flow in D^^t) using the 



vorticity-streamfunction formulation 

dtuj + [i^,uj]^0, (2.1) 
ou = Aijj. (2.2) 

The velocity is given by {u, v) = V"'"'0 := {—dyip, dxi/j), uj — 'V'^-v :— d^v—dyU 



is the vorticity, and [f,g] V / • Vg' = d^f dyQ — d^gdyf is the Jacobian. 
Steady flows satisfy ["0, =0. 

A convenient way of of defining the domain boundary dD^^tj is as the level set 
B{x,y; A{et)) = of some prescribed function B. Since dD^^t) is a material 
curve, 

dtB + [ip, B]^0 for B{x, y; A{et)) = 0. (2.3) 

Assuming that V-B on dD\, this can be inverted to give the boundary 
condition 

ilj = eb{et) on dD^^,t). (2.4) 
Since ip is determined only up to an additive constant, we set 

/ b{et) dl = 0. (2.5) 

It is clear from (2.3)-(2.4) that t/j is proportional to dA/dt on ODa- 
Exploiting the smallness of £, we expand the vorticity and streamf unction in 

u — ou 

Our aim in this section is to compute the leading-order flow ip^'-'\ given its 
initial value and the boundary deformation b{6t), and to show that it depends 
only on the boundary shape A and not on its time history. 

First we note that the boundary conditions (2.4) imply that on dD^, 

jjj^^^ = b and z/'^") = forn = 0, 2, 3, • • • . (2.7) 

Since the total vorticity uj is advcctcd by the flow, the boundary dD^ is a 
vorticity contour and thus we can take cu = u^'^^ and o;*^"^ = for n = 1, 2, • • • 
there. 

Substituting (2.6) into (2.1), we find 

[^(°),a;(°)] + dtJ'^ + e [^(^\a;(°)] + e [^(°\a;(^)] + + 0{e^) = 0. (2.8) 



5 



If the fluid flow is stable in the absence of boundary deformation, we expect 
that the flow will evolve only slowly when the boundary is deforming. We 
therefore introduce the slow time 

T = et, (2.9) 

in terms of which (2.8) becomes 

[^(0)^^(0)] + ^Q^jo) ^ ^ [^(1)^^(0)] + £ [^(0)^,^(1)] + 0{e^) = 0. (2.10) 

At leading order we obtain 

[^(o),a;(o)] = 0. (2.11) 

Taking into account the fact that ^^'^^ = on dD, we flnd that the leading- 
order flow ^^'^^ is instantaneously steady. The relation (2.11) implies that there 
exists a scalar function G relating cu^^^ and ^^^\ 

iP^""^ ^G{uj^'''>;t). (2.12) 

As noted, the function G depends on the slow time r, regarded here as a 
parameter for reasons which will be apparent later. We define F as the inverse 

of G: G{F{u; r); r) = u for every u and r. With an abuse of notation, we will 
often write G' for G'ocuW = W^^VVcj^") and F' for F'oV^W = Vcu(°VW^°^; 
what is meant will be clear from the context. 

At 0{e) we have 

+ [^(^),a;(°)] + [^^^\uj^'^] = 0. (2.13) 
Using (2.12), the second term can be written as 

[^^'\u;^'^ = G"[u;(°), AV'W] = G'AV'W], (2.14) 
with which (2.13) becomes 

9^a;(°^ + [(/.,a;(°)] =0, (2.15) 
= [1-G"A]V'^'^. (2.16) 

These two equations imply that the leading-order vorticity u^^^ is rearranged 
by a divergence-free velocity field V"'"0 with related to the first-order stream- 
function V'^^^ by (2.16). 

2.2 Determination of the Eulerian flow 

We now show how the leading order flow can be determined from (2.12) 
and the fact that Lo'^'^\t) is a rearrangement of its initial value uj^^\Qi). We 
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make the following two assumptions on ip^^")-. 



HI. The leading-order streamfunction tp^^'' is such that it has only one critical 
point in D\ (which is necessarily elliptic) and is nonlinearly stable in the sense 
of Arnold. 



We recall that Arnold stability [cf . iHolm et al.l (119851 )] requires that the steady 



streamfunction ip^^^ satisfies either (i) < Ci < G" < C2 < oo, or (ii) < 
1/cpoi < Ci < —G' < C2 < 00. In the second condition, Cpoi is the Poincare 
constant for the domain D/y, namely the smallest eigenvalue fi of the problem 

(A + n = in Da with u = on dD^. (2.17) 

These conditions ensure that the steady flow is either a minimum or a maxi- 
mum of the energy for fixed vorticity distribution. Note that HI implies that 

-Cpoi<F'<oo, (2.18) 

a condition which will be useful below. The assumption HI is stronger than 
that made in WV but it considerably simplifies the solution of (2.40)-(2.42) 
below. 



For the second assumption, we need a little more notation. Let s denote a 
variable conjugate to ijj^^^ in Da, satisfying [?/^'^°\ s] = 1. Denoting the differ- 
ential arclength along the curve ip^^^ = const by dl, we have ds = d//|V?/'°|. 
We then assume: 



H2. There exists a > such that, for all values of c assumed by ■ip^^\ 

ds < — . (2.19) 



In the context adiabatic invariance, this condition is natural: the left-hand 
side of (2.19) gives the period of rotation of fluid parcels along the streamline 
= c; its boundedness ensures that a time-scale separation between this 
period and the time scale of the boundary deformation exists for sufficiently 
small e. As noted in WV, H2 holds if cua 7^ at the fixed point of ip\. 

For concreteness, we henceforth assume that, at t = 0, our domain is param- 
eterised by Aq and we choose our coordinates in C such that Aq = 0. Fur- 
thermore, we fix in Dao = -Dq a steady leading-order flow 'ip^^\x,0) = ipoi^) 
satisfying H1-H2. Considering only the leading-order flow ip^^\x,t) for the 
moment, we then claim that, assuming H1-H2: 
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PI. The How i/j^^^ is uniquely determined by (i) the shape of the deformed 
domain Z)a, (ii) the steadiness condition (2.12), and (Hi) the fact that the 
vorticity cu^'^^ = Aip^^^ is obtained by rearrangement of the initial vorticity 

UJQ = Al/jQ. 

As a result, the leading-order flow at a fixed time t depends only on the shape 

of the deformed domain at t (parameterized by A{et)), and not on the history 
of shapes at intermediate times (parameterized by the path A(r), < r < 
et). One may draw an analogy with adiabatic invariance in finite-dimensional 
Hamiltonian systems with slowly-varying parameters: here the amplitude is 
completely determined by the instantaneous value of the parameter, not by 
its time history. 

To emphasize the fact that the leading-order fiow depends on A instanta- 
neously, we introduce the notation 

^(0) ^ and cu(°) = ua (2.20) 

for the leading-order streamfunction and vorticity. We also write G{-;t) —: 
Ga{t){-) and F{-;t) =: Fa(,-)(-). These define the scalar functions G\ and Fa, 
both of which have A as a parameter. When (and only when) no confusion 
may arise, we will often write G\ o cua as G\, G'^ o u!\ as G'^, and similarly for 
Fa. 

Like the other results in the present paper, the claim PI is only local: it holds 
only for sufficiently small domain deformations. A similar result is proved in 
WV with a precise functional setting and a different (weaker) set of hypothe- 
ses. The main idea, which we repeat here for reference, is to reformulate the 
problem in terms of the area-preserving diffeomorphism 

qa: Dq ^ Da: X ^ qaX, (2.21) 

which effects the vorticity rearrangement, that is, such that cua — ^^o^ gj}- In 
terms of the pull-back 

g*A : C{Da) C{Do) : / ^ / o ^a, (2.22) 
this can be rewritten as 

ijA = (^a')*^o. (2.23) 

Note that, for fixed t^o and ua, is not defined uniquely by (2.23): rearrange- 
ments along the lines of constant vorticity clearly have no effect. Correspond- 
ingly, the time derivative of qa is not necessarily the divergence-free velocity 
field V^cf) appearing in (2.15)-(2.16), but the equality 

^qaX = V^[(t){gAX] A) + w{ujA{gKx))\ (2.24) 
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holds, where vj is an arbitrary function of one variable. This non-uniqueness, 
of no importance as far as t/'A and cl^a are concerned, will play a crucial role 
when particle positions are examined in §3. 

The map gp^ satisfies a nonlinear partial differential equation obtained as fol- 
lows. Since a;A is a steady flow in Da, we have using (2.12), 

a;A = A(G'Aoa;A), (2.25) 

so applying we find 

a;o = (7lA(^A')*(GAoa;o). (2.26) 

The associated boundary conditions are 

gA(dDo)^dDj, and G^^Gq (2.27) 

for the boundary value of uq, the latter following from the fact that — 
ipQ = and u!\ = ujq on dD\. 

The partial differential equation (2.26), with g/^ and Ga as unknowns, is shown 
in WV to have a locally unique solution (modulo translations along vorticity 
contours) using a contraction mapping argument. This establishes PI and pro- 
vides a way of computing g\ and G\, and hence u!\ and ipA. Alternatively, PI 
can be established using the stability assumption HI: the associated charac- 
terisation of steady flows as energy extrema makes it clear that the steady flow 
V'A is the (locally unique) extremum in with vorticity distribution fixed by 

(jjQ. 

We now consider the infinitesimal version of (2.26), that is, we consider the 
change in g\ corresponding to an infinitesimal deformation of the domain. This 
yields a different construction for g/>^, based on integration over A, and provides 
all the ingredients needed for the computation of the first-order correction to 
■0A and of the Lagrangian flow. 

2.3 Infinitesimal deformations 

In what follows, we will often make use of the fact that many important 
quantities are linear in the boundary deformation rate A := dA/dr. We will 
regard these as resulting from the pairing between the vector A e T\C and 
a differential one-form belonging to a dual space. For instance, the function b 
appearing in the boundary condition (2.4) is linear in A. This makes it possible 
to define the one-form 

/3(-;A):Ta>C^C(9L>a) (2.28) 
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by 

b^f3-A, (2.29) 

where ■ denotes the pairing between vectors and one-forms. Working with 
differential forms of this type gives a compact notation, factoring out the 
dependence in A. At the same time, it allows for a geometric interpretation of 
our results as explained in the Discussion. 

Let d be the exterior derivative in C Since 

-^gA = d<7A ■ A, (2.30) 
and g\ is area preserving, we can define a function- valued one-form $ by 

d(7A = o (2.31) 
An equivalent statement is that 

d^I/ = [*,/] (2.32) 

for any A- independent function /. The initial domain Aq = and initial flow 
flow ipo having been fixed, $ depends only on A and takes its value in the 
space of functions in D\; explicitly, 

$(-;A) :rA£^C(L>A). (2.33) 

Alternatively, since D\ C M^, we can think of $ as a map from to the 
cotangent bundle of C, that is, 

$ : ^ T*C. (2.34) 



From (2.31), $ can be recognized as a connection one-form encoding the 
change in g\ that results from an infinitesimal change in A (making use of 
the identification between functions and divergence-free vector fields on D^; 
see the Discussion for a more precise interpretation) . Introducing (formal) co- 
ordinates {A^} in (a neighbourhood of A = in) £, $ takes the more exphcit 
form $ = $mdAm, where each is a function in D\ (sum over repeated 
indices is implied here and henceforth). 

Taking the exterior derivative of (2.23) gives 

dwA+[*,^A] = 0, (2.35) 

after using the definition (2.31). In components, this reads 



dAm 



+ [$„,u;a]=0. (2.36) 
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Unlike u^, the leading-order streamfunction is not simply rearranged as A 
changes. Applying d to -^a = Ga(<^a) and using (2.35) show that 



dV'A + [$,V'A] ^dGAOcuA. (2.37) 

Here and elsewhere, in dGA o i^a the exterior derivative is taken with respect 
to the (parametric) dependence of Ga on A, so d(GA o'^a) = dG^A oi^a + {G'^ o 
(^k) da;A. 

From (2.35) we derive a property of $ which will be useful later. Taking d 
of (2.35) and using the fact that d^ = 0, leads, after a short computation 
detailed in Appendix A, to 

d^+\[^ ^^^wouK (2.38) 

for some two-form w ouj\: Tj^C C{Da). Here, the bracket [■ A ■] is defined 
in coordinates by 

[a A 7] = [am, 7n] dA^ A dA„ (2.39) 
for any two one-forms a = OmdAm and 7 = 7„dA„ with values in C{Da_). 
(Note that in contrast with [/, /] = for any function /, [a A «] 7^ in 
general.) Equation (2.38) can be interpreted as an integrability condition, an 
infinitesimal version of the statement that g\ is a unique function of A modulo 
displacements along lines of constant u!\. 

We next obtain a dynamical equation for $. Taking the derivative d of (2.26) 
and after a little algebra, we find this equation in the form 

(A - f;) V^a] - A(dGA o u;a) = 0. (2.40) 

The corresponding boundary conditions are 

$ = /3 on ODa, (2.41) 

which follows from (2.4), (2.29) and (2.16), taking into account that 

dGA o cuA = on ODa, (2.42) 

which follows from the fact V'a = on ODa- 

Equation (2.40) is the infinitesimal version of (2.26). It can be solved for [$, iJja] 
and dG^A as follows. For any function u e C{Da), we define the projection Pau 

by 



where, as before, ds = dZ/jV-i/'Al- It follows from this definition that if u is 
constant on a contour of constant ipA, Pau = 0. Note that since the contours 
of uja and ipA coincide, a equivalent definition of Pa could have been given 
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in terms of integrals along vorticity contours ci;a = c. The regularity of Pa is 
guaranteed by H2. Letting (p = [ip/y, <I>] — dG^ ° ^a, and using the facts that 
Pa[iI>A: = [V'A, and PA(dGA o <^a) = 0, we then have 

[V'a,^] = Pa¥' and dGA o a;A = (1 - Pa) V'. (2.44) 
Hence we can write (2.40) as 

(A-F1Pa)(^ = 0, (2.45) 

which is a linear pseudodifferential equation involving (p only. Following (2.41) 
and (2.42), the boundary conditions for cp are 

ip=[iJA,l3] ondDA. (2.46) 

It is shown in Appendix B that (2.45)-(2.46) can be solved uniquely for (p. Us- 
ing (2.44), we recover [ip^, $] and cIGa- From [t/'Aj *^']! $ can be inferred up to 
an arbitrary function- valued one- form depending on x through ^a or, equiv- 
alently, through uja- Thus, there is an equivalence class of forms $ satisfying 
(2.40) which is associated with the gauge transformation 

^ ^ <l> + UouA, (2.47) 

where noa;A is any function-valued one-form depending on x through uj\. This 
non-uniqueness simply reflects at the infinitesimal level the non-uniqueness of 
g\. A particular, uniquely-defined representative of the equivalence class of 
$ could be taken to be Pa^, that is, the one with vanishing average along 
streamlines ipA — const. This is an arbitrary choice, however, and we will see 
in the next section that another choice imposes itself naturally. 

As shown in WV, once we solve the hnear problem (2.40), the solution of 
the nonlinear problem (2.26) for qa follows, at least in a neighbourhood of 
A = and subject to sufficient smoothness of the flow and the domain. For a 
fixed sequence of boimdary deformation, that is, for a fixed path A(r) C C, 
one can in principle solve (2.40) and find $ for each A as long as the fiows 
■ipA encountered along the path satisfy H1-H2. Note that (2.40) is consistent 
with the integr ability condition (2.38) which thus remains satisfied along the 
path (see Appendix A). This confirms our main conclusion, namely that qa is 
independent of the path, up to translation along contours of cua- 



3 Eulerian Flow: First-order Correction 

We now turn to the derivation of the first correction ijj^^^ to the leading-order 
fiow ipA- This derivation is necessary, in particular, to determine the trajecto- 
ries of fluid particles over the 0{e~^) time scales of interest. Remarkably, t/j^^^ 
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can be derived from the knowledge of $ alone. In the process, the gauge of $ 
is fixed in what we argue is a natural manner. 

The derivation starts by noting that (2.24) and (2.31) imply that and $ • A 
differ by a function of a;A only. Thus we can write 

= $*-A, (3.1) 

where 

= + n* o u;a- (3.2) 

Here $* is a specific member of the class of equivalent one-forms it corre- 
sponds to the unique choice of the gauge 11 = 11* in (2.47) that ensures that 
(3.1) holds. In this sense can be seen as a natural choice of connection form. 
The computation which follows shows how it can be obtained from Pa$. 

Since like (f), is linear in A, we can write 

^(1) ^ ^(1) . (3 3) 

where, \E'*^^\ like is a function- valued form; their relationship follows from 
(2.16) as 

il-G'^A)¥^^ = Pa'^ + U'' oujy. (3.4) 

Both ^'^^^ and 11* can be deduced from (3.4). To show this, we make use of 
the constraint imposed by the material conservation of the total vorticity cu. 
This imphes that the function 

A{fl;u}) :— / d^x — area bounded by {cu — fl}, (3.5) 

where int{c<j = Q} denotes the interior of the curve = f2, is an exact invariant 
of the dynamics. Since both the total and leading-order vorticities a; and u!\ 
are rearrangements of the initial vorticity cuq, 

A{n; Lu) = A{n; uja) = A{n; uJo). (3.6) 

Expanding u — ua + sw^^^ + • • ■ , we then have 

A{n; UA + eJ^^) - A{n; ua) - 0{e^), (3.7) 

or, after some manipulations, 

/ u^^^ ds = 0. (3.8) 

The last equation, which can be rephrased as (1 — Pa) uj^^^ = or 

(1 - Pa)A*W = 0, (3.9) 
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provides the solvability condition for (3.4). Indeed, applying (1 — Pa) to (3.4) 
gives 

U* oujA = {1-Pa)¥^\ (3.10) 

This reduces (3.4) to 

(A - f;Pa)*« = -F;Pa$. (3.11) 
The associated boundary condition follows from (2.7) as 

= /3 on ODa. (3.12) 



Equation (3.11) is well posed, with a right-hand side that is uniquely defined 
in spite of the gauge freedom in The operator on the left-hand side is the 
same as that in (2.45) and hence its invcrtibility can be established using the 
same arguments, detailed in Appendix B. Once is determined from (3.11), 
n* follows from (3.10), and the natural connection $* is obtained. Note that 
(3.2) and (3.10) imply that it satisfies 

(1 - Pa)*^^^ = (1 - Pa)$* ^ / ds = / ds. (3.13) 

This relation turns out to be crucial for the computation of fluid particle 
trajectories in the next section. 



4 Lagrangian Flow: Geometric Angle 



In this section we study the evolution of fluid (or tracer) particles in our flow 
over a timescale r = 0(1). 



4-1 Hamiltonian Formulation 



Up to this point, our description of the Eulerian dynamics has been (mostly) 
coordinate- independent. But in order to study particle positions, we need to 
introduce explicit coordinates (,t, y) in D\; {x, y) is chosen to coincide with 
the fixed coordinates in the ambient space through which D^iet) moves. 

The evolution of a particle with position {x{t),y{t)) moving with the fluid is 
governed by the Hamiltonian system 

dx di^ dy dijj 
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with the streamfunction ijj acting as the Hamiltonian. Our aim here is to 
obtain an estimate of {x{t),y{t)) with an error of 0{e) for r = 0(1), so in the 
rest of this section we put 

H{x, y, t) = ^A(.t)(a:, y) + e^'^^\x, y; A{et)) ■ A (4.2) 

in place of in (4.1), keeping in mind the validity of this approximation. 

For e = and hence A constant, the Hamiltonian (4.2) is integrable. For e 0, 
two types of perturbations make it non-integrable: the slow time dependence 
of ip\ introduced by the time dependence of A, and the 0{e) change introduced 
by the addition of -ip^^^ = ^l/^^-* ■ A. We examine th e comb i ned e ffect of these 
two perturbations following closely the approach of iBerrvl (jl985h . 

Since the leading-order Hamiltoniar i iJj\ is integr able for fixed A, we first 
change to action-angle variables [cf. (lArnoldl . Il989l . pp. 297ff)]. At each {x,y), 
we define the action / by 

H^,y) = ^ / , ^ ^dxdy=: ^A(^a). (4.3) 

The angle 9 is defined as the variable conjugate to J, [1,0] = 1. It is 27r- 
periodic since the contours of ip^ are closed and it is related to the variable s 
used earlier by 

27ids = A'{ijA)de. (4.4) 

The canonical transformation (x, y) ^ (J, 6) is obtained by a generating func- 
tion S{I, y; A), with 

X = -— and u = -—. (4.5) 
oy ol 

Solving these implicit equations, we can write 

x = X{I,e;A) and y = Y{I,e;A). (4.6) 
With these and (4.3), we define 



MI; A) = V^(°)(X(/, 6; A), ¥{1, 6; A); A) = A-\27rI; A), 
9; A) = ¥'\XiI, 9; A), F(/, 9; A); A). 



(4.7) 



Here and in the rest of this section, we denote by a hat quantities considered 
as functions of (/, 9). 

So far we considered a fixed value of A. Now let A evolve slowly in time, 
A = A{et). The equations of motion in (J, 9) variables are 

dJ dH , d^ dH ,^ , 
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where the new Hamiltonian H{1 , 9; A) is related to H{x, y; A) by 



H{i, e- A) = H{x{i, e- A), e- A); a) + ^. (4.9) 

We note an abuse of notation here: properly speaking H = H{x,y,t), but 
since the t-dependence only enters through A{et) and its derivative, we have 
written H = H {x , y; A{et)) . Differentiating the definition 

S{I,9;A) = SiI,Y{I,e;Ay,A) (4.10) 

with respect to t at fixed (/, 6) gives 

dA dS dSdY dS ^^^^ dA 

dS = — + - — = +XdY (4.11) 

dt dt dy dt dt dt ^ ' 

Upon substituting dS/dt into (4.9), we obtain that 



H{I,e;A) = Ml)+^'^'^'\l,0;A)-A 

+ e{ds{i, e- A) - x{i, e- A)dr(/, ^; A)} ■ A. 



(4.12) 



Since particles are attached to contours of vorticity u = const, which only 
deviate by 0{6) from the corresponding contours of cja, the action can only 
vary by 0{e) over timescales r ~ 0(1). This is also evident from direct com- 
putation: since H is independent of 9 at leading order and is periodic in 9 at 
the next order, 

^ = -.^{^«+d5-XdF}.A, (4.13) 



and the principle of averaging [cf. (lArnoldl . Il989l . § 52)] implies that / changes 
only by 0{e) for r = 0(1). 

The behaviour of the angle variable is more interesting. From (4.8) we have 

The change in the angle := 9{t) — 9{0) can then be expressed as = 
A6'dyn + A6'geo- The dynamic phase A6'dyn simply arises from the instantaneous 
frequency of the particle, which is the first term in (4.14) above, 

A^dyn = -QJ '^Hr')il) ^t' . (4.15) 

The other terms make up the geometric angle A^geo- 
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4-2 Geometric angle A^, 



In thi s subs ection we sho w that, as in the finite-dimensional cases of iHannay 
fll985h and iBerrvl dlQSSh . the angle A^geo can be understood in geometric 
terms as the (an)holonomy of a connection as a closed path is traversed in a 
parameter space. 

From (4.14), the geometric angle A^geo can be written as 
A^geo = r tttI^^'^ (/, e; A(r')) + dS{I, 6; A(r')) 

(4.16) 

- XiI,e;Air'))dYiI,e;Air'))} . — dr' 



This form suggests that A6'geo depends only on the path traversed in C and 
not on its time parametrisation. The terms inside the braces do depend on / 
and 6, but as shown earlier, the total variation of the action / is of 0{e) over 
the timescale of interest. The dependence on the periodic variable 6 can be 
removed by averaging. For any function / periodic in 6, let 



1 



(4.17) 



We note that by (4.4) this is essentially equivalent to the projection 1 — Pa. 
Applying (■) to (4.16) and replacing I{t) by /(O), we find 



^rnlc ^^^'^ + dS-XdY) + Oie) 



(4.18) 



where C\ is the path traversed in C and where the integrand depends only on 
/ and A. 



Because of the arbitrariness in the angle coordinates (depending on our choice 
of ^ = for each A), the geometric angle is only unambiguously defin ed when 
the p ath C \ is cl o sed, t hat is, when A(r) = A(0) = 0. Following iHannay 
( 119851 ) and iBerryl (jl985l ). we consider this scenario, which is illustrated in 
Figure 2. Since is exact, it vanishes when integrated around Ca- Using 
Stokes' theorem in C, the remaining terms in (4.18) can be written as 



— I ld¥^^ - dx A dr 

Ol JSa ^ 



(4.19) 



where is a two-dimensional surface bounded by C\. The second term is 
identical to that obtained by Berry (1985; eq. (18)) for general Hamiltonian 
systems depending slowly on time; the first term results from the 0{e) change 
to the Hamiltonian induced by the boundary deformation. 
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Fig. 2. Angle change for a cyclic domain deformation. As A(r) describes the closed 
loop Ca C with interior Pa, the fluid domain D\ is deformed and returns to 
its original shape. Fluid particles remain on vorticity contours which approximately 
coincide with streamlines. The position of the particles along vorticity contours is 
defined by the angle-like variable 9 whose total change includes the geometric 
contribution A^geo which depends only on the geometrical properties of Pa- 

Now (^(^)) = ($*) by (3.13), so the first term in the integral can be written 
as 

(d^W) = (d$*). (4.20) 
As for the second term, we use the fact that is a canonical transformation 
(since it is area-preserving) to write 

X(7,^;A)=^aX(7,^;0), (4.21) 

thus defining the transformation to action-angle coordinates for all values of 
A in terms of the transformation at A = 0. 



It follows from this and (2.31) that 

dX(7,e;A) 

We then have 



dX dY 

dXAdy = — — — - dA^ AdA„ = 

dAm OAn 

= if$ $ 1 dA A dA 



X(/,e;A) 



(4.22) 



dA™, A dA, 



dY dX 
1[$A$] = A<l>* 



(4.23) 



where the last equality follows from the fact that the bracket is independent 
of the gauge choice for $. Furthermore, the fact that the transformation to 
action-angle variables is canonical implies that 



A$*] = A**], 



(4.24) 



where the second bracket is in terms of {1,0). We can therefore write (4.19) 
as 



(4.25) 
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A last step, detailed in Appendix A, shows that 

K := d<l>* ^ •l'^] = d^'^ + A (4.26) 

This function-valued two-form can be recognised as the curvature of the con- 
nection and according to (2.38), it depends on space through uj\ or, equiv- 
alently, through / only. This property stems from the constraints of area and 
vorticity preservation imposed on the particle motion along vorticity contours. 
The average in (4.25) is therefore superfluous, and we obtain the result: 



P2. The geometric angle of a particle caused by the slow deformation of the 
boundary is given by 

A^geo = 4? f (4.27) 

where k, given in (4.26), is the curvature of and depends only on the action 
I of the particle and on the domain shape parameterised by A. 



5 Small boundary deformation 



As a concrete illustration of the developments so far, we consider the case 
where the total boundary deformation is small. In general, computing g\ in 
(2.21) requires either solving (2.26) or integrating the differential equation 
(2.15) with boundary condition = 6 [cf. (2.41)], which would have to be done 
numerically. Analytic progress is possible, however, if one considers boundary 
deformations that are sufficiently small for a perturbative approach to be 
applicable. I n this section we develop suc h an a pproach systematically using 



Lie series [cf. iLichtenberg and LiebermanI (119921 )]. 



Let 6 he a formal small parameter characterising the smallness of the boundary 
deformation. The function B{x] A) defining D\ can then be expanded as 

B{x; A) = Bo{x) + 6B,{x- A) + d^B^ix; A) + ■ ■ • , (5.1) 

where Bq is independent of A (recall that Aq = 0), Bi is linear in A, B2 
quadratic, etc. Since is close to the identity and area-preserving, it may be 
regarded as the fiow at 'time' S of an associated ^-dependent divergence-free 
vector field given by V"'"p for some function p{x, 6). (Here p is to be considered 
to live in C{ND x /), where ND C is large enough to contain all relevant 
-Da and / C M.) Correspondingly, the pull-back of defined in (2.22) satisfies 
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Expanding p in powers of 5 as 

p = pi + 5p2 H (5.3) 

and introducing into (5.2) lead to the expansions 

gl^l + •] + y ([p2, •] + [Pi, [Pi, •]]) + • • • , (5.4) 

{9K r = 1 - •] - y ([P2, ■] - [Pi, [Pi, •]]) + • • • . (5.5) 

Computing g*p^f for an arbitrary A-independent / using (5.4), taking the ex- 
terior derivative and identifying with (2.32) lead to 

$ = 5dpi + -(dp2 + [dpi, Pi]) + • • • , (5.6) 

up to an arbitrary function of a;A, after using the Jacobi identity. 

Introducing (5.4)-(5.5) into (2.26) leads to a sequence of partial differential 
equations for the coefficients of p. The ffist two read 

(A-F^)[pi,V'o]-Axi=0, (5.7) 
(A - Fi)[p2,M - 2Ax2 = -2[pi, A[pi,^o]] + [pi, [pi,^o]] (5.8) 

+ A[pi, [pi, V'o]] + 2[pi, Axi] - 2A[pi, xi], 

where Fq is shorthand for Fq o ^q, and we have introduced the expansion 

Ga o u;o = -00 + ^Xi + 5^X2 > • • • , (5.9) 

with Xi-i ^ — 1)2- ••, depending on x through ujq{x). These equations are 
supplemented by the boundary conditions 

[pi,Bo] = -Bi, (5.10) 
[p2,i?o] = -2i?2 - 2[pi,5i] - [pi, [pi,i?o]], (5.11) 

to be applied on the curve Bq{x) = B{x;0) = 0. The formulation is then 
relatively simple, with all the equations to be solved in the original domain 
Dq. The functions Xn^ ^ = 1) 2, • • • are found from solvability conditions. These 
can be made exphcit using the same method as in the treatment of (2.40). For 
instance, using the projection operator Pq associated with lines of constant 
ipo, (5.7) can be rewritten as 

(A-F^Po)v^ = (5.12) 
where (f :— [pi, -^o] — Xi? implying that 

[pi,M = and xi = (1 - Po)^- (5.13) 
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Once the p„, n = 1,2, are computed, the leading-order vorticity and 
streamfunction foUow readily from 

OJA = Uo- S[pi,UJo] - -j{[p2,uJo] - [pi, [pi,'^o]]) H , (5.14) 

V'A='0O-(^([Pl,'0o]-Xl) 

- y([P2,^o] - [pi,[pi,^o]] -2x2) +••• ■ (5.15) 



To find the first-order correction to the Eulerian fiow, (3.11) needs to be 
solved by expansion in powers of 5. This is conveniently done by pulling back 
this equation to the original domain Dq. To do this, we define the pull-backs 
(denoted by overbars) and their expansions as 

fr(i) := = + 52^(1) + . . . ^ (5.16) 
|) := g,* $ = S^i + 6^^2 H 

= 5dpi + -(dp2 - [dpi, Pi]) + ■■■ , (5.17) 

where the last equality follows from (5.4) and (5.6). Introducing these pull- 
backs into (3.11) and noting that 

12VX1 



= F^(^o) - ^[^o(^o)]^^ + 



leads to 



(A-F^Po)^S') = -F^Po«>i, (5.18) 
(A - F^Po)^'^'^ = -F^Po5>2 + A[pi, - [pi, AfrW] 

-F'^a¥^\ (5.19) 

These equations, involving the same invertible operator as (5.12), can be solved 
to find with ^^^^^ deduced after application of (gX^)*- The natural gauge 
of $ then follows from (3.2) and (3.10). Alternatively, one can first compute 
the pull-back !»*, which is obtained from the relations 

= Po$ + n* o a;o and H* o t,;o = (1 - Po)^^^^ (5.20) 
inferred from (3.2) and (3.10), and then deduce by pushing forward with 

To compute the curvature k and the geometric angle, there is in fact no need 
to push forward ^^^^ and indeed, from (4.7) and (4.21), we see that ^^^^ 
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and '^^^^ are related by the A-independent transformation 

^(^)(/,^;A) = ¥^\X{I,e;0)]A) (5.21) 

defining the action-angle variables in the original domain Dq. Since and 
obey an analogous relation, they are essentially equivalent: in particular, 
dl>* = dl>* and [l>* A $*] = [$* A $*]. The curvature k, in (4.26) can therefore 
be computed directly from i>* in a straightforward manner as 

«; = d$*- A$*]. (5.22) 

Note that, in principle, the first two terms in the expansion of or need to 
be computed in order to obtain a leading-order approximation to the geometric 
angle. This is because <l>^ is independent of A, d$* = and hence n = 0{6^). 
The computation can however be shortened by observing that the average of 
l>2 along streamlines, that is, (1 — Po)^2' o^^Y 0(5'^) quantity genuinely 

needed if the averaged form (4.25) of k, is used. In turn, (1 — Po)^2 can be 
approximated by (1 — Po)^'2^\ as the pull-back of (3.11) indicates. The latter 
quantity satisfies a relatively simple equation, obtained by applying (1 — Pq) 
to (5.19) to find 

(1 - Po)A^'^^) = (1 - Po){A[pi,^«] - [pi, A'I'W]}, (5.23) 
after using (1 — Po)A^^i^ — which follows from (3.9) at leading-order in 5. 



6 Nearly axisymmetric flow 

We now consider a simple example where the computations of g\ and other 
relevant quantities can be carried out explicitly to 0(5^). We assume that 
for A = 0, the fiuid domain is the disc (r, a) e [0, 1] x [0, 27r]. The deformed 
domain is defined by 

r = 1 + 5 ^ A^e""'^ - ^ ^ |A„|2 + 0{5'), (6.1) 

TO TO 

where the A^ G C satisfy A*^ = A_„j, with * denoting complex conjugate. 
The multi-dimensional parameter A is therefore infinite dimensional: A = 
{A^ G C : m G Z}. Area preservation at 0(5^) requires that Aq = and the 
introduction of the 0(5^), cr-independent terms. 
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6.1 Arbitrary axisymmetric flow 



The initial flow is taken to be axisymmetric, with vorticity 

uJo{r) = ^(r^o(r)) , 

where the prime denotes differentiation with respect to r. For this flow, (5.7) 
reduces to 

^Ad^p, + 2(^^)'(dlp^ - Up,) + l{rx[)' = 0, (6.2) 

with Xi a function of r only. The corresponding boundary condition is obtained 
from (5.10) in the form 

'9.Pi = -E^-e''"'' atr = l. (6.3) 

m 

The solvability condition for (6.2), found by integration with respect to a G 
[0,27r], imposes that (rx'i)' = 0; boundedness of Xi then implies that Xi is a 
constant which we can take equal to zero: Xi = 0- 

The solution of (6.2) for pi is then found as the Fourier series 

pi(r,a) = EA„Pl,„(r)e""^ (6.4) 

m 

with p*„ = pi-m- Equation (6.2) does not constrain the m = mode pi,o; 
this is the result of the gauge freedom for gp^. A convenient choice is 

Pi,o = 0. (6.5) 
Introducing (6.4) into (6.2) gives a second-order equation for pi^m, namely 

, , / ,, 1 , 2 - \ / , 1 \ , , 

y^o[Pl,m - -Pl,m + P^'^ ) + (^Pl,m " -pl,m j = 0, (6.6) 

with associated boundary condition 

Pi,m = — at r = 1. (6.7) 



m 



There is a close connection between this equation ar id the Rayleigh equatio n 



for the normal modes of axisymmetric flows (e.g., iDrazin and Reidl . Il98ll ): 
(6.6) can be recast as the Rayleigh equation for zero-frequency modes, with 
r~^ifjQPi^m as the unknown function. Of course the non-homogeneous boundary 
condition for pi ^ differs from the homogeneous boundary condition usually 
considered for the Rayleigh equation. The connection is useful nevertheless: 
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the absence of zero-frequency normal modes that can be estabhshcd from the 
Rayleigh equation when ipQ ^ (as guaranteed by the hypothesis H2) imphes 
the existence of a unique solution to (6.6). 

We note that the solution for the m — 1 mode, which describes a rigid trans- 
lation of the disc, is independent of ipQ and given by pi,±i = lAj-ir. Not 
surprisingly, this corresponds to a uniform displacement field V"'"pi. 

The vanishing of indicates that the vorticity-streamf unction relationship 
is unchanged at leading order in S. This is a particularity of axisymmetric 
flows which makes it worthwhile to carry out the calculation to 0{6^) so as to 
demonstrate how a non-zero X2 obtained; this is described in Appendix C. 

With pi determined by its Fourier series (6.4), l»i is given by 

TO 

Because pi = and 1 — Pq is simply the average along circles, (5.18) indicates 
that (1 - Po)^P = 0. Equations (5.20) then imply that $^ = $1. In other 
words, our choice (6.5) provides the leading-order connection with its natural 
choice of gauge which corresponds to vanishing average along the circles r = 
const. Expanding ^^^^ in Fourier series as 

*?^ = E^S!mWe'"^'^dA^' (6-9) 

TO 

(5.18) is reduced to the set of ordinary differential equations 

-Uj'M}n=-^'oPl,m (6.10) 

with = 0. The associated boundary conditions are found from (3.12) as 

= - at r = 1. (6.11) 

m 

Solving (6.10) gives the first-order correction "^^i^ to the Eulerian fiow to 
leading order in S. 

As discussed at the end of §5, the computation of the geometric angle to 
leading order requires not only $^ but also $2 '^^^ minimise computations, 
(1 — Po)'^'2^- This is deduced from (5.23) which reduces to the ordinary dif- 
ferential equation 

lii^ii^ - Po)*^^^} = (1 - Po){A[p^,m - [Pi, A*«]}. (6.12) 
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Solving this equation leads to an expression for (1 — Po)^2 ■ Taking the dif- 
ferential yields the first term of the curvature k in (5.22) as 

(d$*) = (1 - Po)d¥^^ + 0{5^). (6.13) 

Note that since pi is linear in and ^i^'* is A-independent, (1 — Po)^2^'* is 
linear in A„j; furthermore, because the averaging 1 — Pq along circles eliminates 
all products in the right-hand side of (6.12) except for those of complex- 
conjugate Fourier modes, (1 — Po)^'2^'* is a linear combination of terms of the 
type A^dAjJ^. Therefore, (d$*) is given by a A-independent linear combination 
of the two-forms dA^ A dAj;^. 

The second term in (5.22) is also 0(5^) and is readily computed from (6.8). 
Averaging along circles gives it the same form as that of (d<l*). This leads to 
the geometric angle in the form 

A^geo = 5' E fm{r)An + 0(5^), (6.14) 

m>0 

for some functions fmir)- Here we have defined 

Am = -^ [ dA^AdA;; 

2 Jva 

which be recognised as (minus) the oriented area enclosed by the path de- 
scribed by A^ in the complex plane. (A positive Am is associated with a 
rotation of the fluid domain in the positive sense.) Unsurprisingly, at leading 
order, the geometric angle is the sum of separate contributions of each Fourier 
mode of the boundary deformation. 

6.2 An example: flow with power-law radial dependence 

As a simple example of an axisymmetric flow, consider the streamfunction 

Mr) = with < q; < 2, (6.15) 

for which (6.6)-(6.7) can be solved explicitly, leading to 

^l,m = Pl,m = , (6.16) 

m 

where 

Pm = ttm — a + 2 and = y/vn? -I- — 2a. 
The leading-order vorticity in the deformed domain is then found to be 

a;A(r,(7) =u;o(r) -(5u;^,(r)EA„.r^--V"''^ + 0(0, (6.17) 
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Fig. 3. Vorticity wa (grey scale) and streamfunction V'A (white lines) of the lead- 
ing-order (steady) flow in a slightly deformed disc. The top left panel m = shows 
the undcformcd flow, with tpQ = r^/"^-^ the other panels show the flows obtained when 
deforming the disc by a single Fourier mode m according to (6.1) with 5|A^| = 0.05. 

with a similar expression for t/'a. Figure 3 shows these approximations to uja 
and t/'a in domains deformed by a single Fourier mode m, with m ranging 
from 2 to 6, in the case a = 1/2. 

Equation (6.10) can also be solved explicitly, with the result 



*i!m = - Tmr^"* + (1 - 7m)r"'" where 7^ = (6.18) 



1 



m 



a 
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Introducing into (6.12) gives 

+ (l-7jF(a„,/?Jr'^-+^'"-2 

where we have defined 



(6.19) 



(6.20) 



Prom this and (6.13) we deduce the first component of namely 



m>0 



+ (1 - 7m)i^(am,/3™)r-"'"+^--'J dA^ A dA;;, + 0{5^). 

(6.21) 

Using (6.16), the second component of k, is found directly to be 



A2 /.27r _ 

([!•* A = ^ / [$1 A $i] d(7 + 0(5^) 

ZTT Jo 



m>0 



(6.22) 



Combining these results with (4.25), and noting that the action-angle variables 
in the undeformcd domain are simply (7, 9) — (r^/2, a), leads to the geometric 
angle in the form (6.14) with 



fm{r) = - (pmr'^--' + qmr"--^^--') , 



(6.23) 



where 



Pm = ^mE{Pm, Pm) - 2/?m(/?m - 1) and Qm = {I - '^m)E{am-, Pm)- (6.24) 

Figure 4 shows the functions fm{r) for m = 2, 3, • • • , 6 in the case a — 1/2. 

A spot check for our results is provided by the limit a ^ 2, corresponding to 
a flow with a uniform vorticity AA. Assuming that = if m 7^ ±2 and 
that A-i-2 trace the unit circle in the complex plane, the domain deformation 
is simply the rotation of a small-eccentricity ellipse with semi-axes 1-1-25 and 
1 — 2S. There is an exact analytic solution for such a uniform-vorticity flow 
in a rotating ellipse; both the direct use of this solution and (6.14) yield the 
same r-independent value for the leading-order geometric angle that appears 
for one full rotation of the ellipse, namely 167r(5^. See Appendix D for details. 
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Fig. 4. Functions fm{^) giving the geometric angle in (6.14) for the axisymmetric 
flow with streamfunction '0o = r^^"^ in a disc deformed by a Fourier mode m. 

7 Discussion 

In this paper, we have used differential geometry mostly as a notational tool, 
its main application being to make explicit the linear dependence of various 
quantities on A. It is nonetheless clear that the objects we are dealing with can 
be given an interpretation in a more abstract geometric setting and that the 
problem may be placed in the framework of geometric mechanics. This would 
be a valuable undertaking, but one which is probably quite involved due to 
the infinite-dimensional nature of the problem, and which is certainly beyond 
the scope of the present paper. However, we now discuss briefly and informally 
the geometric context of our results in order to elucidate the meanings of the 
connection one-form $*, the form of its curvature k, etc. 

Central to our development is the group G = SDiff(M^) of area- preserving 
diffeomorphisms of the plane. Taking the initial domain Dq as a reference 
domain, the subgroup H <Z G which maps Dq to itself, viz., 

H^{geG:gDo^ Do}, (7.1) 

is of particular importance. In terms of G and H, our parameter space C, 
the space of all possible shapes A of the domain Da, can be realised as the 
quotient G/H. Indeed, each right coset of in G contains all diffeomorphisms 
mapping Dq to D^ for a particular A, since any two such diffeomorphisms g 
and g' are related by g' = gh for some h E H. We can therefore identify L 
with G/H. Another important subgroup is 

Hq ^ {h e H : ujqo h ^ ujq} (7.2) 

containing all area-preserving diffeomorphisms in Dq which leave the initial 
vorticity distribution invariant. A rearrangement ujqo g oi the initial vorticity 
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ujq can be identified with an element of G/Hq. 

Armed witli tliis setup, we can interpret our results geometrically. The key 
point is to regard G and G/Hq as principal bundles, both with G j H L 
as base manifold. Finding the (leading-order) Eulerian flow for each domain 
shape A then amounts to flnding a lift from GjH to G/Hq\ flnding the La- 
grangian particle position amounts to flnding a lift from G/ H to G\}\ Propo- 
sition PI, stating that the leading-order Eulerian flow depends only on the 
domain shape, says that the lift from G/H to G/Hq is path- independent; in 
other words, it deflnes a section of GZ-f^oH] In contrast, the lift from G/H 
to G, which gives an approximation to the particle position, depends on the 
path in G/H and in fact on the speed with which the path is traced. There 
is, however, a contribution that is independent of speed; for cyclic domain 
deformations, it is quantifled by the geometric angle given in P2. 

It is worth commenting on the meaning of the one-form that appears in the 
geometric angle. One way of deflning a lift in a principal bundle is by means of 
a vector- valued one- form, i.e. a linear map from T{G/H) to TG, describing the 
vertical (along-flbre) displacement associated with any given displacement on 
the base manifold G/H. Such a form can be recognised as a connection form. 
In our context, TG is the space of divergence-free vector flelds over M^, which 
can be identifled through the use of a streamfunction with the space of real- 
valued functions C(M^). Thus, a lift can be deflned by a connection one-form 
over G/H C with values in C(]R^). This is precisely the interpretation we 
give to $*. With the geometric interpretation of $*, the subsequent results are 
clear: (4.26) is the standard expression for the curvature of $*, the geometric 
angle (4.27) is given by the holonomy of and the standard conclusion about 
the geometrical angle in flnite-dimensional systems is recovered. 

Throughout the paper, we emphasize that our results are local in nature: the 
prediction of an instantaneously steady Eulerian flow, for instance, holds only 
if the domain deformation is such that HI and H2 are always satisfled. This is 
necessarily the case if they are satisfled initially and the domain deformation 
is sufficiently small, but it may well continue to hold for larger deformations. 
It is nonetheless interesting to speculate about the dynamics when either HI 
or H2 fails in the course of the evolution. If HI fails, the flow ceases to be 
Arnold stable and likely becomes spectrally unstable. We can then expect the 
flow to become highly unsteady and, in the absence of dissipative mechanisms, 
remain so regardless of subsequent deformations of the domain. 

^ In our geometric description the interior and exterior of D\ are treated on the 
same footing; this can be done because our formulas, with suitable boundary con- 
ditions as I a; I ^ oo, would also apply to a fluid flow outside Da- 
^ We stress again the local nature of PI: globally, there are many possible steady 
flows for a given vorticity distribution and a given domain D\. Geometrically, this 
implies that the section of G/Hq is multivalued. 
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The failure of H2, on the other hand, corresponds to the appearance of stream- 
hnes for which the orbiting period of particles becomes large. When the pe- 
riod becomes comparable with the time scale of the domain deformation, our 
asymptotic approach clearly breaks down. This can happen when the flow is 
driven by the domain deformation towards a change in topology, with the cre- 
ation of hyperbolic stagnation points and separatrices. How the flow evolves in 
this situation is unclear, but some understanding could be gained by investi- 
gating the problem where the initial steady flow ipo h as a hyperbo l ic stag nation 
point (and satisfies HI — the Kelvin-Stuart vortex in lHolm et al.l ( 119851 ) is just 
one example). This problem can be viewed as a g eneralisation of th e classical 
critica l-layer problem for parallel shear flows (e.g.. lStewartsonl (119781 ) : lMaslowe 
( 1l986l ) and references therein). In this generalisation, the separatrix plays the 
role of the zero- velocity critical line (along which H2 is obviously violated); 
by analogy, it can be expected to be also surrounded by a narrow critical 
layer where complicated nonlinear dynamics occurs. We plan to investigate 
this problem in future work. 



We conclude this paper by remarking on the possible extension of our results 
to flows in three dimensions. In three dimensions, the dynamics of an inviscid 
and incompressible fluid is determined, as in two dimensions, by a form of con- 
servation of vortic ity, although in this case i t is as a vector that the vorticity is 
transported (e.g.. [Arnold and Khesinl . Il998l ). This suggests that our approach 
for the determination of the leading-order Eulerian flow in deforming domains 
can be adapted to the three-dimensional setting. The technical conditions for 
the well-posedness of the equations for g/<^ are however likely to be significantly 
more complicated than in two dimensions. 



The evolution of the fluid-particle position seems, at first sight, to pose a 
very different problem in three than in two dimensions, since the velocity 
field is divergence-free and not Hamiltonian. However, particle trajectories for 
(non-Beltrami) steady solu tions of the Euler equations are known to be inte- 
grable dArnoldl . Il965l . Il966f ) because they are confined to surfaces of constant 
Bernoulli function (Lamb surfaces). There is, therefore, a simple characteri- 
sation of the fluid-particle positions in steady flows, analogous to the action- 
angle characterisation in two dimensions. This could be used for slowly time- 
dependent flows to quantify the effects of a cyclic boundary deformation as 
was done to obtain the geometric angle in this paper. 



A general difficulty with three-dimensional flows, however, is the absence 
of gener al stability results similar to those obtained by the energy-Casimir 
method (IHolm et al.l . Il985l ). Instabilities cannot therefore be excluded (on the 
contrary, they are the rule rather than the exception), and the effect of their 
competition with the slow evolution of the leading-order flow would need to 
be assessed carefully. 
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A Properties of $ 

In this Appendix, we give details of the derivations of three useful expressions. 

1. We first derive (2.38). Applying d to (2.35) gives 

= d^a;A = -d[$, loa] 

= -[d^,ujA\ + A da;A] (A.l) 
= -[d$,a;A]-[$A[$,u;A]], 

where the bracket [• A •] is defined in (2.39). Now 

A [^,UJa]\ = -[^rn, [^n,i^A]] dA^ A dA„ 

= [ua, + [^^A, [^m, ^n]]} dA^ A dA„ (A.2) 

= [$A[$,a;A]] + [u;A,[$A$]] 

(we have used Jacobi's identity, [/, [g, h]] + [g, [h, /]] + [h, [/,(?]] = for any 
three functions /, g and h, to arrive at the second equality); therefore 

-[$A[$,a;A]] = |K,[$A$]] (A.3) 

and from (A.l) we find 

[d$+i[$A$],u;A] =0 

^ d$-F A$] = wocua (A.4) 

for an arbitrary function- valued two-form w o uja- 

2. Next we show that (A.4) can be derived directly from (2.40). This provides 
a consistency check for our developments. Let us first compute 

d[$,^A] = d{[^m:ijA\ dA^} 

= [d$^, ijA] A dA^ + d^A] A dA^ (A.5) 
= [d$, Va] - a dV'A] 
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and 



d(dG'A o ua) = -{dG'j, o wa) A duA = {dC^ o ua) A [$, lua] , . „s 

where the last equahty can be verified by computation in coordinates. 
Writing (2.40) as 

A [$, ^a] - [$, ooa] - A(dGA o a;^) = (A.7) 

and taking d, we find 

= A d[$, V'a] - d[<^, cua] - Ad(dG'A o ua) 
= A[d$, iPa] - A[<l> A dtpA] - [d$, uja] + [$ A dc^A] + A[$ A d^A] 
= (A - f;) [d$, V^a] + a [$ a (dG-A - d^A)] + [$ A dcuA] 
= (A - FX) [d$ + |[$ A $], ^a]. (A.8) 

A couple of identities have been used to arrive at the last equation. The first 
one is 

[$A[$,V'A]] = i[[$A$],V'A], (A.9) 
which is proved in the same way as (A. 3). The second identity is 

A da;A] = A [$, Fa o ^^]] = [$ A F;[$, ^a]] . , . 

= F;[$A[$,V'A]]-[$,n]A[$,V'A]=n[$A[$,V'A]], ^ ■ ' 

where wc have used [$, F;^ o ■i/^a] A [$, V'a] = i^A V'a] A [$, V^a] = for the 
last equality. 

The desired result (A.4) is recovered by noting that the operator (A — F^) is 
invertible by hypothesis and that [d$ + |[$ A $], ■?/'a] = on ODa- The latter 
can be verified by difi'erentiating (2.35) and evaluating it on ODa- 

3. Finally, we establish the formula 

d$-|[l>Al>] =d$ + |[$A$]. (A.ll) 

Its application to the natural connection $* shows that k, is independent of 9. 
Our proof starts by noticing that [$ A$] = [$ A$] because the transformation 
to action-angle variables is canonical. Thus (A.ll) is equivalent to 

d$ = d$+[$A$]. (A.12) 
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This is established by direct computation as follows 

d$(/, 9] A) = d$(X(J, e- A); A) 

= d^{x- A) + ^ (dX ■ V$„) dA„ 

n 

= d$(a;;A) + 5^[$^,$„] 



(a=,A) 

dA™ A dA, 



, , s m 

(a3,A) 



d$(x;A) + [$ A$] . (A. 13) 

(a!,A) 



B Solution of (A -F'^P^)u = f 



Here we show that the problem 



(a-f;pa)7^ = / 

r] = g on dDx 



(B.l) 



has a unique solution 77 when F;^ > — Cpoi everywhere in as follows from 
the hypothesis HI. 

We start with an identity. Let u and v be such that Pam = and Pa^" = v. 
We have 

j uvdxdy — J ^(j) uv ds^ djpA — J u i^ji v ds^ d^pA — 0. (B.2) 

From this it follows that the projection Pa is orthogonal in L'^{D\), in the 
sense that for any (sufficiently smooth) function w 

dxdy = ^{|PA«^r + 2(Pa«^)[(1 - Pa)H + 1(1 - Pa)^'!'} dxdy 
^ l^{\PAw\^ + \{l-PA)w\^}dxdy. 

(B.3) 

Using (B.2), we find that the operator (A — -Fa Pa) is self-adjoint for any 
functions u and v which vanish on dD, 

[ v{A- Fj^PA)udxdy = [ u{A- Fj^PA)v dxdy. (B.4) 
Jd Jd 

Moreover, (A — F'^^ Pa) is coercive under the hypothesis F' > — Cpoi- To show 
this, we first combine (B.3) and Poincare inequality to obtain 

[ \Pau\^ dxdy < [ \u\'^dxdy<—f \Vuf dxdy (B.5) 
Jd Jd Cnn\ Jd 
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for any function u vanishing on dD. It is then clear that 

J^u{A-Fj^ Pa) udxdy = - ^{l + F^iPAU^} dxdy<0 (B.6) 
when Fj^ > Cpoi everywhere in D, with equahty obtaining only when u = 0. 

Returning to the problem (B.l), we extend g to cIDa and let fj = i] — g. The 
problem thus becomes 

(A - F;^Pa) fi = -Ag + F'j^P^g (B.7) 

with boundary conditions r/ = on dD\. We have shown that the operator 
(A — -FaPa) on the left-hand side is self-adjoint and its associated bilinear 
for m is coercive, so by the Lax-M ilgram lemma (assuming compactness, etc., 



cf. iGilbarg and Trudingeii (119771 )) a unique solution rj exists for (B.l). 



C Second-order terms in nearly cixisymmetric flows 



At order 0(5^), p2 is found from (5.8) to satisfy 



(C.l) 

= 2[pi, A[pi, V^o]] - [Pi, [Pi, AV^o]] - A[pi, [pi, i)^]]. 
The boundary condition (5.11) can be written as 

dap2 = d„ (drpid^pi) - {d^pif + Yl at r = 1, (C.2) 

m 

after some manipulations. The interest of this form is that, when (6.7) is 
taken into account, it is clearly consistent, with both sides having a vanishing 
cr-average. A solvability condition for (C.l) is obtained by averaging over a, 
leading to 



\27cr 



(draPlf + -^{dlaPlf - -dapApi 



da] . 



This equation determines X2 uniquely up to an irrelevant arbitrary constant. 
When it is satisfied, (C.l) can be solved for p2, yielding a solution in the form 
of a Fourier series 



m 
m 



with P2m — P2-m and p2,o = 0. The functions p2,m satisfy an inhomogeneous 
version of (6.6) obtained from (C.2); clearly, they are quadratic in the A^. 
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D Rotating ellipse 



Consider a fluid inside an ellipse with semi-axes a and b that is rotating around 
its centre with a (possibly time-dependent) angular velocity eA(r). The equa- 
tion of the ellipse is given by 



5(x,r) = ^ + |-l 



0, 



where 



Xi = X cos A + ?/ sin A, 
X2 = —X sin A + ?/ cos A 



are Cartesian coordinates in a frame rotating with angular velocity eX. An 
exact solution for the fluid motion in such a rotating elhpse is provided by the 
uniform- vorticity flow with streamfunction 



a? 



+ 



2{p? + 62) 



(x\ 



Xr 



where K is a constant. In this streamfunction, which can be verifled directly 
to satisfy the boundary condition (2.4), the flrst term can be identifled with 
the second with and there are no higher-order terms in e (see, e.g., 
Jeffreys and Jeffreys! . Il974j . p. 421, for a derivation). Action-angle coordinates 
(/, 9) for this flow satisfy 



Xi 
X2 



2Ia/b cos 9, 



'2Ib/a sine. 

In terms of these variables, the streamfunction, or Hamiltonian, becomes 



H{I, 



2KI 

ah 



+ 



eXI \a? 



ah 



+ 62 



[a cos 9 — h sin 



[a^ cos^ 9 + h'^ sin' 



(D.l) 

In this expression, the last term between round brackets comes from the time- 
dependence in the (canonical) transformation from {x,y) to {1,9); it simply 
corresponds to a rigid-body rotation with angular velocity — eA. The geomet- 
ric angle is derived by writing the evolution equation for 9, averaging, then 
integrating in time. This gives 



geo 



AA 

2ab 



fa2 - b 



2\2 



- {a^ + h' 



a^ + h^ 

where AA is the total angle rotated by the ellipse. 



(D.2) 
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Note that the averaging is in fact unnecessary, since the Hamiltonian H does 
not depend on ^ as a simpUfication of (D.l) indicates. We do not perform 
this simpUfication here, however, in order to retain the two terms in (D.2) 
separately. This facihtates the comparison with the general formalism of §4. 
The first term in (D.2) stems from the correction in the streamfunction tp^^^ 
(which here corresponds to a potential fiow), while the second stems from the 
slow time-dependence of the leading-order fiow; thus these two terms can be 
identified with the two contributions (d$*) and — A in (4.25). 

We can use the exact formula (D.2) to verify the approximate results for 
slightly deformed axisymmctric fiows obtained in §6. A rotating ellipse of small 
eccentricity is represented by the deformed disc (6.1) with A-|-2(''") describing 
a unit circle in the complex plane and all the other equal to zero. The 
corresponding semi-axes are then 

a = 1 + 25 + 0(0 and 6=1-25 + 0(0, 

with the 0{6^) corrections ensuring that ab = 1. Introducing this into (D.2) 
and considering a full rotation AA — 27r gives the geometric angle 

325' - (2 + 165')] + 0(5^) = -2tt + 16nS^ + 0(5^). (D.3) 

An equivalent result is obtained from the developments in §6. Since the only 
non-zero parameter A^ is A-1-2, (6.14) reduces to 

A^geo = 5V2(r)A + 0(53). 

The uniform- vorticity corresponds to the limit a — > 2 in (6.15), so that a2 — 
P2 — 2, 72 = 1/2, and /2(r) = 8, independent of r. Because a full rotation of 
the ellipse is obtained when A2(t) covers twice the unit circle, A2 — Stt and 
hence 

with contributions 32ttS^ and —16tt6^ from (d<l>*) and —i([<l>*, <!>*]), respec- 
tively. The discrepancy of — 27r when compared with (D.3) results from a 
different definition of the angle 6, which is measured from an axis rotating 
with the ellipse in the calculation leading to (D.3) while it is measured from 
a fixed axis in §6. 



A^geo = TT 
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